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Abstract. - We present an experimental study of the outflow of a hopper continuously vibrated 
by a piezoelectric device. Outpouring of grains can be achieved for apertures much below the 
usual jamming limit observed for non vibrated hoppers. Granular flow persists down to the 
physical limit of one grain diameter, a limit reached for a finite vibration amplitude. For the 
smaller orifices, we observe an intermittent regime characterized by alternated periods of flow 
and blockage. Vibrations do not significantly modify the flow rates both in the continuous and 
the intermittent regime. The analysis of the statistical features of the flowing regime shows that 
the flow time significantly increases with the vibration amplitude. However, at low vibration 
amplitude and small orifice sizes, the jamming time distribution displays an anomalous statistics. 
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Introduction. — Since the end of the Middle Ages 
sand clocks have been customarily used to measure time 
flowing (see for example an historical monography in [1]). 
Nowadays, hoppers are widely used in the industry to store 
and deliver granular matter [2] . In both situations one uses 
a striking property of granular matter as the flowing rate 
does not depend on the height of material in the container 
but rather, on the size of the orifice. This was rationalized 
in the 60's by Beverloo, who proposed a simple scaling law 
for the discharge rate essentially based on the free falling 
velocity at a scale comparable to the aperture [3,4]. The 
scaling law includes an empirical modification to account 
for small orifices (sizes of the order of a few grain diame- 
ters). In this case, one observes blocking of the flow due 
to the presence of static vaults forming at the exit [5]. 
This can be problematic for industrial processes since it 
severely limits the possibility to deliver small quantities of 
material in a reliable way. In recent years this issue has 
been studied in detail with the aim to characterize more 
thoroughly the jamming probability and the flow prop- 
erties near this limit [5-9]. One way to deal in practice 
with the problem of blockages is to add vibrations in order 
to restore the flow by breaking the arches [10]. Previous 
works have been devoted to the influence of vibrations in 
a hopper discharge [11-15]. In most cases, however, the 
main interest was the influence of vibrations on the flow 
rate. As the motivation was usually a practical situation, 
these studies often deal with strong vibrations and large 



orifices, well above the jamming limit. Here we focus on 
small orifices, which are prone to jamming. Indeed, this 
question is not just a technical one, but relates to the fun- 
damental properties of granular matter and can be put in 
framework of current interrogations on the jamming tran- 
sition of soft matter [16]. In the present case, deep inter- 
rogations arise about the capability for a granular packing 
to flow in a restrained geometry and on the true nature of 
the block/unblock mechanisms. From a theoretical point 
of view this set of questions can be framed in the con- 
text of the rigidity transition driven by the existence of a 
marginal mechanical state , i.e., where the average num- 
ber of contacts is just sufficient to ensure a solid structure 
to the packing (the isostatic limit [17]). In this case, one 
observes a large density of so-called '"soft modes'" corre- 
sponding to collective granular motion at low energy cost. 
At the outlet of a hopper the confining pressure vanishes, 
thus giving a situation where, if of any importance, anoma- 
lous critical behavior should be revealed (at least from a 
finite size scaling perspective). In this paper, we present a 
novel insight on the problem, with an experimental design 
of a hopper that allows a fine statistical study of granular 
flows in this limit. The continuous variation of two crucial 
parameters (diameter opening and vibration amplitude) 
allows to deal with the probability of arch formation and 
thus, to go well below the usual jamming limit. In this 
way, we are able to investigate the properties of the gran- 
ular flow down to the physical limit of one grain aperture. 
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Experimental setup. — The silo is made of a plexi- 
glass cylindrical tube (inner diameter 1.9 cm) and is specif- 
ically designed to control continuously two important pa- 
rameters that will allow an extensive study of the transi- 
tion to a blocking state, namely, the aperture opening and 
the vibration amplitude. To this purpose, the design of the 
bottom outlet is special. The cylinder's bottom edge was 
obtained by a planar cut at 45 dcg. The resulting ellip- 
tic area is partially closed by a flat rectangular surface on 
which a thin piezoelectric transducer is mounted, its di- 
ameter being larger than the silo diameter. The blocking 
piece can be moved along the cut in the direction of the 
larger diameter of the ellipse via a micrometric translation 
stage (see Fig. [TJ. The hole hence produced, is character- 
ized by a minimal opening distance L (see inset of Fig.Q]). 
The piezoelectric transducer is driven with a sinusoidal 
vibration at a frequency of 350 Hz and it vibrates in the 
direction perpendicular to its surface. Note that the vi- 
bration is directly applied to the granular material inside 
the silo, but not to the wall, which is not touching the 
vibrating transducer. The granular material consists of 
glass beads, their diameter being 1.55 ± 0.05 mm. Prior 
to the experiments, we measured the energy input due to 
the vibrations induced by the piezoelectric device. An ac- 
celerometer was buried at the top of the fully filled silo and 
for different electric voltage amplitudes, the r.m.s. accel- 
eration was monitored, with the outlet closed. A linear re- 
lation is obtained between the electric voltage and the ac- 
celeration in the bulk of the granular column. We checked 
that the acceleration docs not change significantly when 
the outlet is open. The typical measured accelerations are 
only a fraction of the gravity acceleration. For the maxi- 
mum induced vibration amplitude (a rms acceleration of 
about 1 m/s 2 ), the typical kinetic energy is much smaller 
(by about five orders of magnitude) than the potential en- 
ergy at the scale of a grain diameter. The mass flux is 
monitored by recording the mass of grains that fall from 
the silo on an electronic scale interfaced with a computer. 
To have a good determination of the time lapses when the 
orifice is jammed or when it flows, a plastic sheet with a 
microphone attached to it is mounted below the silo out- 
let. Thus the sound produced by the impact of the falling 
grains can be monitored and time lapses inferred. 

Dependence of the mean flow with the orifice 
size. — The determination of the mean mass flow Q was 
carried out by measuring the time taken to deliver a large 
quantity of grains (about 2 x 10 4 ) . The conventional pic- 
ture for the outpouring of grains through a hole is embod- 
ied in the so-called Beverloo's scaling [3,4], that relates 
the mean flow rate to the aperture diameter. In this sce- 
nario, grains are released from an imaginary vault having 
the same radius as the exit orifice. In our case, it is not 
clear what is the right length scale A controlling the size 
of the dynamical vault that actually releases the grains 
(and therefore determines the mean flow). So in order 
to benchmark our experimental setup, we measured the 




Fig. 1: Experimental set-up: M, microphone; S, scales; P, 
piezoelectric transducer; T, translation stage. The top of the 
granular column is open at the atmospheric pressure. A sketch 
of the bottom outlet is also shown with the variables L and D. 
The piezoelectric transducer vibrates in the direction perpen- 
dicular to its surface as indicated by the black arrows. Note 
that the walls of the hopper are not touching the piezoelectric, 
but only the grains. The top right insets show the signal as 
measured from the microphone (top, signal level in arbitrary 
units) and the binarized signal (bottom). 

mean flow Q (in the absence of vibration) as a function of 
the size of the exit orifice for the cases where no perma- 
nent jams are observed. We have taken as a characteristic 
length the projection of L (see Fig. [I) on the horizontal 
axis: D = ^L. To check whether the Beverloo scaling 
is valid, the length scale A = (Q<7 _1 ^ 2 p -1 ) 2 ^ 5 , where g is 
the gravity acceleration and p the density of the beads, is 
computed. As seen in Fig. [21 the data without any vibra- 
tion are compatible with the Beverloo scaling Q cx D 5 / 2 . 
Remark that we do not intend to discuss here the valid- 
ity of this scaling, which is a priori not obvious for small 
orifices (and especially with the present geometry). This 
has been done elsewhere [18]. The relevant point is that 
although the shape of the outlet is not circular and its area 
does not increase as D 2 , the data follow conspicuously well 
Beverloo's scaling down to the blockage limit. This seems 
to indicate that there is only one structural length scale 
controlling the granular flux, and that it is directly linked 
to the aperture distance D. 

The influence of vibrations. — The applied vi- 
bration is quantified through the rescaled acceleration F, 
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Fig. 2: (Color online) Length scale A normalized to the bead 
diameter d as a function of rescaled aperture size D for various 
vibration amplitudes. The dashed line is a fit of the data for big 
D; the slope is close to one. The arrows mark the minimum 
rescaled orifice sizG Dmin (definition explained in the text). 
Inset: Value of D m in / f d as a function of F. 



which is the rms acceleration normalized by the gravity 
acceleration g. Note that the pertinent parameter could 
also be the kinetic energy. We do not intend to discuss 
here the influence of the driving frequency, and only the 
driving amplitude at a given frequency is changed. 

In the absence of vibration and for small apertures, 
there are situations where the flow is continuous for some 
time and then it can eventually get blocked forever. The 
application of a vibration can unjam the flow (it depends 
on how energetic the vibration is). With applied vibration, 
the flow can also get blocked, but the outpouring can re- 
sume after a finite time. In other words, an intermittent 
regime is observed. In order to extract rational flow prop- 
erties we must define a strict and precise working protocol 
that involves the finite mass of grains we seek to deliver 
and the time we are ready to wait before deciding that 
the flow is halted. The protocol chosen involves a fixed 
amount of grains (approximately 2 x 10 4 beads); if this 
mass is delivered without any jam that lasts for 100 s or 
more, then the hopper is considered to be flowing. If there 
is a jam that lasts for 100 s or more, the run is finished 
and the silo is considered to be jammed. The time interval 
of 100 s was so chosen because our experiments indicate 
that it is very unlikely that the flow will resume by itself 
without vibration after that time lapse. For a given T and 
D, the mean flow Q is calculated as the mass discharged 
divided by the total time of measurement. The size of the 
orifice D is then decreased and the measurements are car- 
ried out again, as long as the waiting time does not exceed 
the predefined limit. In this way, for each vibration am- 
plitude we obtain D m i n (T), the smallest size of the orifice 
for which the flow has not been halted during a time lapse 
longer than 100 s before the mentioned amount of material 
is discharged (i. e. the smallest value of D for which the 



hopper is considered to be in a flowing state, as obtained 
with the mentioned protocol). In some way, D rn i n might 
be related to the size of the orifice above which jamming 
is extremely improbable in a non vibrated silo [7]. 

The data arc displayed on Fig. [21 where the values of 
D m in for different T are indicated by arrows. The first 
salient feature is that for orifices of sizes above D m i n (T = 
0) (corresponding to the black arrow) the influence of vi- 
bration is negligible, i.e. the mass flux does not change 
whatever the value of T. Second, the vibrations allow the 
possibility of flow below this orifice size. The more ener- 
getic the vibration is, the lower is the minimal value of D 
for which the silo can flow without getting jammed. The 
last striking point is that for a finite value of T, one can 
reach by linear extrapolation the value D m i n = 1 below 
which hard grains cannot go through the aperture (see in- 
set of Fig. [2]). This result proves conceptually the design 
value of such an activated outlet that is able to deliver 
granular flow down to the ultimate physical limit of one 
grain size (D m i n = 1). However, we need now to push 
the study further and characterize the flowing dynamics 
because below the value of D m i n (T = 0) the mean flow as 
defined before does not follow anymore a simple Beverloo's 
scaling. 

Intermittency and Jamming. — In the presence of 
vibrations at the outlet, therefore, the usual limitation of a 
minimum size to deliver grains from a hopper can be over- 
come. However, in this regime the outlet flow displays a 
strong intermittent dynamics that needs to be character- 
ized in more detail. Using the mass recorded on the scale, 
we examine directly the exit flow continuously over time. 
On Fig.[3]four time series corresponding to the mass deliv- 
ered on the scales are represented for four different vibra- 
tion amplitudes, the exit orifice being the same -a certain 
value below D min (T = 0). As the vibration decreases, 
the curves clearly display jam periods of larger duration. 
These events are rare or nonexistent for apertures larger 
than the value of D min (T = 0). 

From the data series as presented in the inset of Fig. [1] 
(impacts of grains on the microphone) and in Fig. [3] (out- 
flow mass delivered on the scale), we propose that this 
dynamical issue can naturally be separated into two sets 
of questions pertaining to the flowing and to the jammed 
regimes: (i)when the hopper is flowing, how long is it going 
to flow and what is the flow rate? (ii)when it gets jammed, 
how long does the blockage persist and what fraction of 
the overall experimental time does it spend in that state? 

These questions are going to be addressed in detail next. 

The flowing regime. — The time intervals Tp during 
which the hopper is flowing are obtained from the impact 
of the out-flowing grains on the microphone (see Fig. [TJ 
inset) for an experiment at constant D and T. As seen 
in the Fig. 2J the distribution of Tp shows an exponen- 
tial tail. Therefore, a characteristic time t(T, D) can be 
defined. The inset shows that this characteristic time in- 
creases with acceleration T and goes to a finite value when 
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Fig. 3: {Color online) The masses delivered by the silo with an 
exit orifice D = 2.42d as a function of time, for four different 
vibrations as indicated in the legend. The dashed horizontal 
lines mark the upper and lower bounds of the zone where the 
mean flow has been measured for comparison (see explanation 
in the text). 



r tends to zero. The fact that one can define a characteris- 
tic value for the flowing intervals Tp is consistent with the 
existence of mean avalanche sizes obtained in [7, 8] in the 
case of a non-vibrated silo where the discharge through 
the small orifice is restored by jet of pressurized air. The 
remarkable point here is that this is also the case when 
vibrations are applied. This hints to an interesting fea- 
ture that will be checked next, namely, that when the 
vibration is applied the restored flow has the same statis- 
tical features that the corresponding flow one would have 
without vibration. As noted earlier, this is true above 
D mi „(r = 0). Let us take a closer look on Fig. [31 which 
corresponds to a diameter D inside the intermittent do- 
main, and in particular, let us consider four stretches in 
the plots where the orifice is not blocked (as marked with 
horizontal dashes). The mean flow for these intervals are 
all similar: 4.05, 3.99, 3.88 and 4.10 ± 0.01 g/s , respec- 
tively (the resolution is the least mean square fit error). 
To investigate this point in more detail, we computed the 
histograms of the flow rate q corresponding to the same 
orifice size for five different vibration amplitudes (Fig. [5]) . 
Note that in this case the flow q is measured as the num- 
ber of beads delivered in one second. The most remarkable 
feature of Fig. O is that the q PDFs mainly consists of a 
Gaussian centered around the same mean value of the flow 
rate whatever T is. The differences appear for the low flow 
rate events corresponding to the small flow rates. It seems 
that a higher vibration amplitude is able to reduce signif- 
icantly those events that seems to bring the system in the 
vicinity of blockage. Consequently, the time of flow can 
be significantly increased as noted before. Interestingly, 
this distribution is in agreement with previous results for 
nonvibrated silos [19]. 
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Fig. 4: (Color online) Histogram of the rescaled time inter- 
vals during which the orifice D = 1.78d has been flowing (in 
semilogarithmic scale). Inset: Mean flowing time as a function 
of the rescaled rms acceleration V. 



The jammed regime. — Let us now explore the jam- 
ming dynamics due to arch formation at the outlet of the 
hopper. These arches are constantly destroyed by vibra- 
tions and then recreated by jamming. The time intervals 
during which the outpouring has been halted are denoted 
by the variable Tj (see inset of Fig. [Q). The PDFs of Tj 
for a fixed D at several T are displayed in Fig. [6] 

A typical time scale to compare with can be obtained 
by calculating the time it takes for a bead to fall its 
own diameter under the acceleration of gravity, namely 
T g = ^2d/g w 0.02 s. For times larger than r g the distri- 
bution decays as a power law with an exponent that de- 
pends on the vibration amplitude (Fig. [6]). Importantly, 
for small holes and low T the power law displays an expo- 
nent smaller than 2 (absolute value), and this means that 
the first moment of the distribution does not converge. In 
other words, the data show that the distribution of jam 
durations can be extremely large and that the mean time 
during which the hopper is jammed will depend on the to- 
tal time that the experiment has been running. Then the 
value of the mean flow can be dominated essentially by the 
longest jam. Due to these features, the term anomalous 
statistics is used here to describe the situation. 

Since the transition to this regime can not be rigorously 
denned in terms of an average flow, we propose instead 
an alternative characterization. The fraction of time $ 
for which the hopper is jammed is computed inside a time 
window of variable length (Tw)- This time window is 
taken inside a long time series during which the experi- 
ment has been running. The measurements are carried 
out for a fixed diameter D = 2.42d and for several values 

of r. 

On Fig[7{a) the values of $ as a function of Tw are 
displayed on a normal-log representation. It can be seen 
that for high acceleration, a steady state is reached after a 
short transient, indicating the existence of a well defined 
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Fig. 5: (Color online) PDFs of the flow rate (obtained by av- 
eraging in a time window of 1 s) with an exit orifice D = 2.42d 
and for five different vibration amplitudes F as indicated in the 
legend. 



Fig. 6: (Color online) Histogram for the time lapses that the 
orifice has remained blocked, in logarithmic scale, correspond- 
ing to D = 1.78d. The dashed line has a slope of two. Note that 
for the smallest acceleration displayed (triangles) the slope is 
smaller than two. 



mean time spent in the jammed state. However, for in- 
termediate accelerations (r = 0.042, 0.029 and 0.017), no 
steady state is reached and the definition of an average 
time for the jams does not make sense. Finally, $ — > 1 for 
r — > because after a finite time the flow is blocked by 
an arch at the outlet, which is not likely to be destroyed 
in the time scales involved in the experiment. Note that 
the slow dynamics (the time is represented in a log scale) 
have some features similar to aging processes. An alterna- 
tive way to describe this transition is shown on Fig. [7] (b), 
where the values of $ at Tyy = 1000 s are represented as 
function of V for D = 2.42d. For the statistical descrip- 
tion, we plot the median (which is a more robust locator in 
the case of anomalous statistics) of 10 runs, and the error 
bars correspond to the median absolute deviation about 
the median, or MAD. As is obvious, the transition is not 
only evidenced by the values of $ going from one to zero 
in the region 0.01 < T < 0.06, but also by the large error 
bars, as expected. 

Now, let us get back to Fig. [7] (a) and in particular 
to the data corresponding to intermediate accelerations 
(r = 0.042, 0.029 and 0.017). In those cases, no steady 
state is reached and 4> presents a positive slope. There- 
fore, it is reasonable to suppose that $ will tend to 1 for 
Tw — > oo. Assuming this supposition, $oo(r) would be a 
step function with two regimes separated by a critical ac- 
celeration r c (see Fig. [7](b)). For V > T c , the mean time 
spent in the jammed state would be well defined, and a 
steady flow rate would be reached. However, for T < T c , 
we could be able to get the silo flowing, but the outflow 
rate would not be well defined, just because the mean time 
spent in the jammed state would not be well defined. The 
existence of a r c is in agreement with the transition to an 
exponent lower than 2 (absolute value) in the power laws 
found for the PDF of Tj (Fig. EJ). 



Conclusions. — We presented an experimental study 
of the discharge of a granular hopper with an external 
vibration applied at the outlet. The design allows a con- 
tinuous control of the vibrational energy injected in the 
granular packing (via a piezoelectric transducer) as well as 
a continuous variation the orifice size. We were thus able 
to measure precisely the size D m i n (T) below which the 
flow is halted and its dependence with the vibration rms 
acceleration T. In the absence of vibration, the grains stop 
flowing below an orifice size D m i n (T = 0). We show that 
-Dmin(r) decreases with an increasing vibration and this, 
down to one grain size, value reached for a finite vibration 
amplitude. This would allow, in principle, a conceptual 
design suited to deliver individual grains from a hopper. 
This surprising feature poses many questions on the fun- 
damental mechanism at work for blockage; at the moment, 
we have no clear theoretical vision or even a scaling argu- 
ment to explain the relationship between D m i n and T. We 
intend to explore this question in future works by changing 
the frequency of vibration, the grain size and the hopper 
geometry. 

Below the D min (T = 0) threshold, we have evidenced 
an intermittent regime. It consists of flowing and jammed 
periods. We describe the flowing intervals as a stochastic 
process for the flow rate. Vibrations only affect proba- 
bilities to observe low flow rates but hardly change the 
distribution around the mean value set only by the hop- 
per aperture. In a future work we will seek to provide a 
quantitative description of the flow statistics and its re- 
lation with the stopping probability. The jammed inter- 
vals are associated with the blockages of the packing, that 
can last a very long time before a continuous flow may 
restart. At low acceleration, we evidence an anomalous 
statistics for the jamming times, characterized by a power 
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Fig. 7: {Color online) The transition from flowing to jammed 
regime, (a) Fraction of the time for which the flow is jammed 
$ as a function of the observation time window Tw ■ Note the 
logarithmic scale on the horizontal axis. Each line is a single 
experimental run. (b) The fraction of time that the system is 
jammed, calculated for Tw = 1000s, as a function of V. Each 
point is the median of 10 different runs, with the error bars cor- 
responding to the median absolute deviation about the median. 
The red dash line represents a hypothetic step function which 
would be obtained for Tw — * oo. The two different regimes 
defined by the step function are named by roman numbers I 
and II. The value given for F c is tentative 



law distribution with no first moment. More precisely, we 
characterized the transition to an anomalous statistics by 
monitoring the fraction of time for which the hopper is 
jammed as a function of the observation time Tw- This 
anomalous statistics is associated with an elusive value for 
the blockage thresholds that may actually depend either 
on the total number of grains available or the time we are 
willing to wait for the flow to resume. This property is 
strongly reminiscent of the anomalous dynamics usually 
observed for creeping flows of glassy materials. It would 
be interesting to relate the flow anomalies described here 
to the new proposition that has emerged recently to de- 
scribe glassy systems and the contribution of the so called 
"soft-modes" dynamics to transport properties. 
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